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ABSTRACT 



By means of an automatic digital computer solutions were ob- 
tained for a number of second order servo problems with backlash in the 
gear train. A complete range of solutions were obtained with variations 
in load friction, motor friction, load inertia, motor inertia, damping co- 
efficient for the combined operation , and backlash angle . These data 
are presented in chart form suitable for the servo design engineer to 
estimate performance of an arbitrary second order system with backlash. 

The program was conducted on the Control Data Corporation Model 
1604 Computer, and this pfogram is presented with explanations of the 
manner in which the differential equations were solved. From data 
generated by the computer phase plane plots are drawn as a further 
illustration of the results of this thesis , and a sample problem is in- 
cluded to show application of the design charts. 

The writer wishes to express his appreciation to Dr. George J. 
Thaler of the Electrical Engineering Department and to Dr. E. J. Stewart 
and Mr. E. N. Ward of the U. S. Naval Postgraduate Computer Center 
for their assistance and encouragement in completing this work. 
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1 . 



Introduction 



In the engineering analysis and application of servomechanisms 
(hereafter called servos) which include gear coupling , the problem of 
backlash must be considered and is often a source of trouble. In theory 
it is desirable to know the system performance of a servo with a given 
distribution of inertia and friction and to accurately determine the size 
of the limit cycle for any given set of variables. In practice, knowledge 
of the performance of a system may be desirable prior to the construction 
of a mock-up or simulation system . 

Existing methods fonthe analytical solution of the generalized 
backlash problem are laborious and complicated and lack of the straight- 
forward approach needed in the preliminary design of a servo . In this 
thesis design charts are presented for the second order servo with back- 
lash to accommodate an arbitrary set of system variables . 

Previous work in this problem has been done at the U. S. Naval 
Postgraduate School by Lutkenhouse in reference a wherein the backlash 
problem was solved by construction of the phase trajectory in the graphical 
phase plane. This thesis extends the work of Lutkenhouse. It is recog- 
nized that the construction of a phase trajectory for a given problem 
requires many hours of tedious calculation and graphical construction. 
Thus, in the collection of a large volume of data required for the con- 
struction of design charts for the second order servo with backlash an 
automatic digital computer was employed . 
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The work of this thesis first required derivation of the system 
equations, then a computer program was developed to adequately des- 
cribe the equations. The burden of the work subsequently fell in pro- 
cessing sufficient data to construct design charts to satisfy the require- 
ments of the problem. 
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2 . 



Equations for the Servo Motion 



in the construction of a set of design charts it is usually desirable 
to provide sufficient flexibility so that the charts may be used for any 
given servo. Thus, equations are presented in Appendix A for the de- 
termination of the system parameters necessary to enter the design charts 
of this thesis for three types of second order servos: the armature control 
motor, the field control motor, and the two phase motor. 

All of these cases behave in a similar manner when coupled to a 
load through a gearbox which has backlash. Initially there is a period 
when the backlash is taken up, and the system performs as a combined 
system. The differential equation describing this motion for an undamped 
natural frequency of one is 

0 + 30 + 0 -1 

where B is equal to twice the damping coefficient, £ 

At such point in the cycle of motion wherein the load velocity 
exceeds that of the motor, backlash starts, and the motor is not directly 
coupled to the load. There are two differential equations describing 
this motion. The equation for the coasting load is 



where D = 



0i. 



JZ 



The equation for the motor is 



where E = J.„ and F = F. , . 

M M 
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This motion in the backlash region continues until the total amount of 
backlash is taken up, and at that point gear contact is made. The action 
of the servo upon contact is not a simple one. This motion can follow one 
of several types of contact where the basic laws governing the motion are 
the law of conservation of momentum, and the law of conservation of 
energy which must both be considered. In this thesis it is assumed that 
the contact is a plastic one with no gear bounce; therefore, there is energy 
loss in the gearbox due to the plastic contact and the law of conservation 
of momentum was applied. Following the gear contact phase the system 
is again combined, and the cycle of events begins again. 

The characteristics of the second order servo with backlash is 
independent of the undamped natural frequency as shown by the develop- 
ment in Appendix A. Therefore, all of the problems worked in this thesis 
and used in the construction of design charts are based on an undamped 
natural frequency of one radian per second, and the step input is equal 
to one radian. Units of radians and seconds are used throughout. 

The amount of backlash used in the program was 0.3 radian, and 
it was a purpose of the thesis to determine the effect of other magnitudes 
of backlash. 
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3. 



Solution of the Equations 



The differential equations for the second order servo with backlash 
were solved on the Control Data Corporation Model 1604 Computer. This 
computer is described in reference b, but briefly, it is an all-transistorized, 
stored program, general-purpose digital computer having a large storage 
capacity and exceedingly fast computation and transfer speeds. Input to 
the computer was by means of magnetic tape and punched paper tape, and 
the output data was recorded on magnetic tape. The data was then printed 
on the International Business Machine 717 high speed printer. 

Figure 1 shows a schematic of the problem. In programming for the 
solution of the equations it was necessary to successively compute a 
point for the combined motion, test this point to determine if the load 
would drift away from the motor, then if so, continue the cycle into the 
backlash region. A second test was then introduced to determine the end 
of the backlash region, and if the backlash was fully taken up the contact 
part of the program was initiated. This complete program was accomplished 
by the program of Table I . 

Several instructions in the program are worthy of mention. First, 
it is noted that a subroutine was used in the actual solution of the dif- 
ferential equations. This subroutine based on the Runge-Kutta Gill numer- 
ical integration is contained in a library tape in the U. S. Naval Post- 
graduate Computer Center, and programming instructions are contained 
in reference c. The program was assembled in accordance with the ”AR" 
assembly routine of references d and e which also contain a full des- 
cription of each instruction used in this program. 
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In the Runge-Kutta Gill numerical integration four successive 
approximations are used in the computation oi a point (in this thesis 
points are CL 01 second apart) . However, it was not desired to print out 
each point, so a count routine was inserted in the program to print each 
eighth point. Thus in the actual printout, an example of which is con- 
tained in Table II, test points are at intervals of .08 seconds. The 
printout format shows the number of the point in column 1 , the time in 
column 1 , the motor speed in column 3, the motor position in column 4, 
the load speed in column 5, and the load position in column 6. The 
number following the figure printed represents a decimal exponent to 
locate the decimal point. 

The remaining item of the program requiring explanation is the 
portion of the program named switch. This routine is inserted to accom- 
modate the various changes in sign caused by the backlash non-linearity, 

A close study of this sub-routine will show that it was not only necessary 
to change the sign of some of the instructions, but it was also necessary 
to change some of the jump instructions caused by multiple sign changes 
in the backlash test. 

With the basic program of Table I stored in the computer, it is 
merely necessary to insert the desired system parameters for the parti- 
cular problem to be solved into the computer locations beginning at 20156 
in order to obtain the solution for any desired servo problem with backlash. 
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4 . Results 

A total of 120 runs for the second order servo with backlash were 
made and the limit cycle lor these problems were determined. Each ran 
required approximately three minutes on the computer and yielded one 
test point suitable for presenting in this report. These test points are 
plotted in Figures 2 through 11, Figures 12 through 15 are phase trajec- 
tories for a number of selected runs plotted from an actual printout. 

Figure 2 is a plot of limit cycle versus friction ratio for a damping 
coefficient of 0.1 and inertia ratios of 1/9, 1/4, 1, 4 and 9. 

Figure 3 is a plot of limit cycle versus friction ratio for the load 
inertia equal to the motor inertia and for damping coefficients of 0.1, 
0.2, 0.3, 0.4, 0.5, and 0,6. 

Figure 4 is a plot of limit cycle versus backlash angle for load 
inertia equal to motor inertia and for a friction ratio of 0.2. This plot 
shows a linear variation for damping coefficients of 0.1, 0.2, 0.3, 0.4, 

0 . 5 , and 0.6. 

Figure 5 is a plot of limit cycle versus backlash angle for load 
inertia equal to motor inertia and for a damping coefficient equal to 0.1, 
This line represents all friction ratios since the - J variation of 
Figure 2 is a non-varying straight line. For a damping coefficient other 
than 0,1 the curve would spread to form a pencil of lines. 

Figure 6 is a plot of limit cycle versus backlash angle for a 
damping coefficient of 0.1 and for a friction ratio of 0.2. Variations 
of the curves of this show inertia ratios of 1/4, 1 and 4. 
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Figures 7 through 11 are plots of limit cycle versus inertia ratio;: 
for friction ratios of 0,1, 0.2, 0.4, 0.6, and 1.0. Each plot shows 
curves for variations of damping coefficient of 0.1 through 0.6. 

Figures 12 through 15 are actual plots of the phase trajectories as 
plotted from the print-out sheets . From these plots can be determined the 
output speed and position with the corresponding time. 

Figure 12 shows a case with friction ratio of 0.2 at an inertia 
ratio of .25 and a damping coefficient of 0.1. Figure 13 is at the same 
friction ratio of 0 . 2 but shows an inertia ratio of 4 and a damping co- 
efficient of 0.6. Figure 14'and 15 illustrates the effect of increasing 
the damping coefficient from 0.1 to 0.2 with load inertia equal to motor 
inertia and at a friction ratio of 1.0. 



8 



5. Discussion 

After the accumulation of the data it was desirable to present the 
information in the most useful and practical form. The size of the limit 
cycle as shown in the figures of this thesis was previously defined as 
deviation from the equilibrium or average value rather than peak to peak 
variation. The data were plotted in various forms in an attempt to find 
the most useful form. Figure 2 is one such presentation where the data 
is plotted for constant damping coefficient versus friction ratio. Figure 
2 was not found to be the optimum form because in subsequent plots at 
different values of damping’coefficient there was appreciable crossing 
and recrossing of lines of constant inertia ratio, and it was found that 
interpolation between these lines was difficult. Further, there was not 
a wide variation of limit cycle with friction ratio as evidenced by the 
level curves of Figure 2. In fact, for the particular case of balanced 
load and motor inertia and damping coefficient there was no variation 
in limit cycle with friction ratio. This plot is shown as a logarithmic 
plot because it was found that the digital program required a number of 
modifications to satisfy the zero friction case. A single case was run with 
zero friction ratio , and the limit cycle was found to be substantially equal 
to that for a friction ratio of 0.04 which is used in this thesis as repre- 
sentative of the zero friction ratio case . 

Figure 3 is a composite plot showing all values of damping co- 
efficient up to 0.6. Figures 4,5, and 6 are presented to show that the 



size of the limit cycle is directly proportional to the amount of backlash. 



Several runs were made with changes in only the backlash to conclusively 
prove this characteristic . Figure 4 shows the variation with damping co 
efficient with all other variables constant. Figures 5 and 6 are similar 
plots with variations in other factors , This linear variation of limit cycle 
with backlash is significant since it means that design charts can be 
drawn for any constant amount of backlash and variation from this value 
is in direct proportion. In this thesis all charts are for a constant amount 
of backlash equal to 0.3 except for those points in Figures 4,5, and 6 
where backlash in a variable . 

Figures 7 through 11 hre the design charts found to be the most 
useful form for design estimations. In these charts the size of the limit 
cycle is plotted versus inertia ratio. The logarithmic scale was chosen 
for two reasons: first, it avoided the zero cases found at the extreme 
ends of the scales and for which a separate program would have been 
required, and secondly, this scale gives a spread to the curves in the 
region of inertia ratio equal to one which appears desirable. 

The design charts presented herein are most useful for the range of 
variables encompassed by the test points. While it is believed that some 
estimates can be made from these charts for systems with very large and 
very small inertia ratios and beyond the limits of the curves shown, best 
results are obtained for values where interpolation between the lines of 
constant damping coefficient is feasible. Appendix B gives a sample 
problem using the design charts . Appendix C gives instructions in the 
use of the design charts . 
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Conclusions 



6 . 

From the results obtained in the accomplishment of this thesis 
it is concluded that the size of the limit cycle of a second order servo 
with backlash can be represented by design charts , one form of which 
are presented in Figures 7 through 11 of this thesis. 

7 . Recommendations 

Inasmuch as it was necessary to complete this thesis during a 
Summer period, time was of essence, and the scope of the program was 
definitely limited. Therefore, the following recommendations are sub- 
mitted as a reasonable extention of the program. 

1 . Rewrite the digital program to include a test to determine the 
actual size of limit cycle for display on the console with no printout 
required. This will permit longer runs with no time required for printout 
of data. Printout capability should still be retained by means of a Jump 
Key, 

2. Extend the program to include the zero friction and zero inertia 
cases , 

3. Extend the program to include a larger range of inertia ratios. 

4. Reprogram the contact phase to include the case of elastic 
contact and to include the "lossy" case for all types of contact between 
elastic contact and plastic contact. 
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9 „ Tables 

TABLE I 

COMPUTER PROGRAM 



REM BACKLASH PROBLEM PLASTIC 

CONTACT RUN 3 
ORG 20000 



20000 


75 4 20022 


START 


SLJ 


4 


INPUT 1 


REMARKS IN THIS 




00 0 00000 




0 


0 


0 


COLUMN 


20001 


75 4 20063 


COMBINED 


SLJ 


4 


RESET 






00 0 00000 




0 


0 


0 




20002 


75 4 60200 




SLJ 


4 


RUNGE 


SET UP GILL 




50 0 00000 




ENI 


0 


0 


(Procedure to solve 


20003 


00 0 20023 




0 


0 


TABLE 1 


differential 




00 0 20034 




0 


0 


DERC 


equations) 


20004 


75 4 60201 


POINTC 


SLJ 


4 


RUNGE +1 






00 0 00000 




0 


0 


0 




20005 


75 4 60200 


LOAD 1 


SLJ 


4 


RUNGE 


SET UP GILL 




50 0 00000 




ENI 


0 


0 




2000G 


00 0 20040 




0 


0 


TABLE L 






00 0 20051 




0 


0 


DERL 




20007 


75 4 60201 


POINTL 


SLJ 


4 


RUNGE +1 






00 0 00000 




0 


0 


0 




20010 


75 4 20130 




SLJ 


4 


COUNT 






50 0 00000 




ENI 


0 


0 




20011 


75 0 20061 




SLJ 


0 


TEST 1 






50 0 00000 




ENI 


0 


0 




20012 


75 4 60200 


MOTOR 


SLJ 


4 


RUNGE 


SET UP GILL 




50 0 00000 




ENI 


0 


0 




20013 


00 0 20023 




0 


0 


TABLE 1 






00 0 20054 




0 


0 


DERM 




20014 


75 4 60201 


POINTM 


SLJ 


4 


RUNGE +1 






00 0 00000 




0 


0 


0 




20015 


75 4 60200 


LOAD 2 


SLJ 


4 


RUNGE 


SET UP GILL 




50 0 00000 




ENI 


0 


0 




2001.G 


00 0 20040 




0 


0 


TABLE L 






00 0 20051 




0 


0 


DERL 




20017 


75 4 60201 


POINT2L 


SLJ 


4 


RUNGE +1 






00 0 00000 




0 


0 


0 




20020 


75 4 20130 




SLJ 


4 


COUNT 






50 0 00000 




ENI 


0 


0 




20021 


75 0 20073 




SLJ 


0 


TEST 2 






00 0 00000 




0 


0 


0 
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20022 

20023 

20024 

20025 

20026 

20027 

20030 

20031 

20032 

20033 

20034 

20035 

20036 

20037 

20040 

20041 

20042 

20043 

20044 

20045 

20046 

20047 



TABLE I ~ Continued 



75 


0 


00000 


INPUT 1 


SLJ 


0 


0 


Provision for changin g 


75 


l 


20022 




SLJ 


0 


INPUT1 


Program if desired 


00 


0 


00000 


TABLE 1 


DEC 




2 


Table for gill combined 


00 


0 


00002 












17 


7 


15075 




DEC 




.01 




34 


1 


21727 












00 


0 


00000 


T 


DEC 




0 




00 


0 


00000 












00 


0 


00000 


UDOT 


DEC 




0 




00 


0 


00000 












00 


0 


00000 


U 


DEC 




0 




00 


0 


00000 












00 


0 


00000 


QU 


DEC 




0 




00 


0 


00000 












00 


0 


00000 


THETADOT 


DEC 




0 




00 


0 


00000 












00 


0 


00000 


THETA 


DEC 




0 




00 


0 


00000 












00 


0 


00000 


QTHETA 


DEC 




0 




00 


0 


00000 












13 


0 


20156 


DERC 


LAC 


0 


B 


Derivative for combined 


32 


0 


20027 




FMU 


0 


U 




31 


0 


20032 




FSB 


0 


THETA 




30 


0 


20164 




FAD 


0 


ONE 




20 


0 


20026 




STA 


0 


UDOT 




12 


0 


20027 




LDA 


0 


U 




20 


0 


20031 




STA 


0 


THETADOT 


75 


0 


60202 




SLJ 


0 


RUNGE+2 


00 


0 


00000 


TABLEL 


DEC 




2 


Table for gill load 


00 


0 


00002 












17 


7 


15075 




DEC 




.01 




34 


1 


21727 












00 


0 


00000 


TL 


DEC 




0 




00 


0 


00000 












00 


0 


00000 


WDOT 


DEC 




0 




00 


0 


00000 












00 


0 


00000 


W 


DEC 




0 




00 


0 


00000 












00 


0 


00000 


QW 


DEC 




0 




00 


0 


00000 












00 


0 


00000 


THETALD 


DEC 




0 




00 


0 


00000 












00 


0 


00000 


THETAL 


DEC 




0 




00 


0 


00000 
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20050 

20051 

20052 

20053 

20054 

20055 

20056 

20057 

20080 

20061 

20062 

20063 

20064 

20065 

20066 

20067 

20070 

20071 

20072 

20073 

20074 

20075 



TABLE I - Continued 



00 0 00000 
00 0 00000 


QTHETAL 


DEC 




0 


13 0 20157 


DERL 


LAC 


0 


D Derivative for load 


32 0 20044 




FMU 


0 


W 


20 0 20043 




STA 


0 


WDOT 


12 0 20044 




LDA 


0 


W 


20 0 20046 




STA 


0 


IEETALB 


75 0 60202 




SLJ 


0 


RUNGE +2 


13 0 20161 


DERM 


LAC 


0 


F Derivative for motor 


32 0 20027 




FMU 


0 


u 


31 0 20047 




FSB 


0 


THETAL 


30 0 20164 




FAD 


0 


ONE 


33 0 20160 




FDV 


0 


E 


20 0 20026 




STA 


0 


UDOT 


12 0 20027 




LDA 


0 


u 


20 0 20031 




STA 


0 


THETADOT 


75 0 60202 




SLJ 


0 


RUNGE +2 


50 0 00000 




EN1 


0 


0 


12 0 20046 


TEST 1 


LDA 


0 


T HEIALB Test to determine 


31 0 20031 




FSB 


0 


THETADOT Start of backlash 


22 2 20012 


JUMP1 


AJP 


2 


MOTOR 


75 0 20001 




SLJ 


0 


COMBINED 


75 0 00000 


RESET 


SLJ 


0 


0 


12 0 20025 




LDA 


0 


T Changes load initial 


20 0 20042 




STA 


0 


TL conditions back to 


12 0 20026 




LDA 


0 


UDOT motor initial conditions 


20 0 20043 




STA 


0 


WDOT if separation has not 


12 0 20027 




LDA 


0 


U occurred 


20 0 20044 




STA 


0 


W 


12 0 20031 




LDA 


0 


THETADOT 


20 0 20046 




STA 


0 


THETALD 


12 0 20032 




LDA 


0 


THETA 


20 0 20047 




STA 


0 


THETAL 


10 0 00000 




ENA 


0 


0 


20 0 20045 




STA 


0 


QW 


20 0 20050 




STA 


0 


QTHETAL 


75 0 20063 




SLJ 


0 


RESET 


50 0 00000 




ENI 


0 


0 


12 0 20047 


TEST2 


LDA 


0 ' 


THETAL Test to determine if 


31 0 20032 




FSB 


0 


THETA backlash has been 


31 0 20166 


BKTEST2 


FSB 


0 


BACKLASH taken up 


50 0 00000 




ENI 


0 


0 


22 2 20137 


JUMP2 


AJP 


2 


CONTACT 


75 0 20012 




SLJ 


0 


MOTOR 



15 



20076 

20077 

20100 

20101 

20102 

20103 

20104 

20105 

20106 

20107 

20110 

20111 

20112 

20113 

20114 

20115 

20116 

20117 

20120 

20121 

20122 



TABLE I - Continued 



36 0 20117 


SWITCH 


SSK 


75 0 20100 




SLJ 


75 0 20107 




SLJ 


50 0 00000 




ENI 


12 0 20062 


LEFT 


LDA 


14 0 20116 




ADD 


20 0 20062 




STA 


12 0 20075 




LDA 


14 0 20116 




ADD 


20 0 20075 




STA 


12 0 20074 




LDA 


15 0 20172 




SUB 


20 0 20074 




STA 


12 0 20117 




LDA 


05 0 00001 




ALS 


20 0 20117 




STA 


75 0 20001 




SLJ 


50 0 00000 




ENI 


12 0 20062 


RIGHT 


LDA 


15 0 20116 




SUB 


20 0 20062 




STA 


12 0 20075 




LDA 


15 0 20116 




SUB 


20 0 20075 




STA 


12 0 20074 




LDA 


14 0 20172 




ADD 


20 0 20074 




STA 


12 0 20117 




LDA 


05 0 00001 




ALS 


20 0 20117 




STA 


75 0 20001 




SLJ 


50 0 00000 




ENI 


00 1 00000 


FACTOR 


OCT 


00 0 00000 






25 2 52525 


ALTERNAT 


OCT 


25 2 52525 






75 0 00000 


PRINT 


SLJ 


12 0 20025 




LDA 


20 0 20252 




STA 


12 0 20027 




LDA 


20 0 20253 




STA 


12 0 20032 




LDA 



0 ALTERNAT Plan for changing 

0 LEFT sign of backlash 

0 RIGHT each cycle 

0 0 

0 JUMP 1 

0 FACTOR 

0 JUMP 1 

0 JUMP 2 

0 FACTOR 

0 JUMP 2 

0 BKTEST2 

0 FAC TOR 2 

0 BKTEST2 

0 ALTERNAT 

0 1 

0 ALTERNAT 

0 COMBINED 

0 0 

0 JUMP1 

0 FACTOR 

0 JUMP1 

0 JUMP2 

0 FACTOR 

0 JUMP2 

0 BKTEST2 

0 FACTOR 2 

0 BKTEST2 

0 ALTERNAT 

0 1 

0 ALTERNAT 

0 COMBINED 

0 0 

0010000000000000 

2525252525252525 

0 0 Print routine 

0 T 

0 BUF 

0 U Motor speed 

0 BUF+1 

0 THETA Motor position 
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20123 

20124 

20125 

20126 

20127 

20130 

20131 

20132 

20133 

20134 

20135 

20136 

20137 

20140 

20141 

20142 

20143 

20144 

20145 

20146 

20147 



TABLE I - Continued 



20 0 20254 




STA 


0 


12 0 20044 




LDA 


0 


20 0 20255 




STA 


0 


12 0 20047 




LDA 


0 


20 0 20256 




STA 


0 


75 4 71000 




SLJ 


4 


01 0 20252 




01 


0 


06 0 00001 




06 


0 


72 0 20126 




RAO 


0 


75 0 20120 




SLJ 


0 


75 0 00000 


COUNT 


SLJ 


0 


12 0 20170 




LDA 


0 


14 0 20171 




ADD 


0 


15 0 20173 




SUB 


0 


22 0 20134 




AJP 


0 


14 0 20173 




ADD 


0 


20 0 20170 




STA 


0 


75 0 20130 




SLJ 


0 


75 4 20120 


OK2PRINT 


SLJ 


4 


50 0 00000 




ENI 


0 


10 0 00000 




ENA 


0 


20 0 20170 




STA 


0 


75 0 20130 




SLJ 


0 


50 0 00000 




ENI 


0 


12 0 20162 


CONTACT 


LDA 


0 


32 0 20027 




FMU 


0 


20 0 20165 




STA 


0 


12 0 20163 




LDA 


0 


32 0 20044 




FMU 


0 


30 0 20165 




FAD 


0 


20 0 20031 




STA 


0 


20 0 20027 




STA 


0 


20 0 20044 




STA 


0 


20 0 20046 




STA 


0 


13 0 20156 




LAC 


0 


32 0 20027 




FMU 


0 


31 0 20032 




FSB 


0 


30 0 20164 




FAD 


0 


20 0 20026 




STA 


0 


13 0 20157 




LAC 


0 


32 0 20044 




FMU 


0 


20 0 20043 




STA 


0 



BUF+2 

W Load speed 

BUF+3 

THETAL Load position 
BUF+4 
DECO 
BUF 
1 

/-I 

PRINT 

0 

INDEX Plan for printing 
INCRONE only every tenth 
TEN point at .08 sec 

OK2 PRINT interval 
TEN 
INDEX 
COUNT 
PRINT 
0 
0 

INDEX 

COUNT 

0 

G Apply law of 

U momentum 

GU 

H 

W 

GU 

THETADOTPuts resulting values 
U of variables in 

W Table 1 and L 

THETALD 
B 
U 

THETA 

ONE 

UDOT 

D 

W 

WDOT 
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20150 

20151 

20152 

20153 

20154 

20155 

20156 

20157 

20160 

20161 

20162 

20163 

20164 

20165 

20166 

20167 

20170 

20171 

20172 

20173 

20174 



TABLE I - Continued 



12 0 20047 




LDA 


0 


THETAL 


20 0 20032 




STA 


0 


THETA 




10 0 00000 




ENA 


0 


0 




20 0 20030 




STA 


0 


QU 




20 0 20033 




STA 


0 


QTHETA 


20 0 20045 




STA 


0 


QW 




20 0 20050 




STA 


0 


QTHETAL 


75 4 20120 




SLJ 


4 


PRINT 




50 0 00000 




ENI 


0 


0 




75 0 20076 




SLJ 


0 


SWITCH 


50 0 00000 




ENI 


0 


0 




50 0 00000 












17 7 66314 


B 


DEC 




.4 


B equals twice the 


63 1 46314 










damping coefficient 


17 7 54000 


D 


DEC 




.125 


D equals the load 


00 0 00000 


* 








time constant 


17 7 56314 


E 


DEC 




.2 


E equals motor inertia 


63 1 46314 












17 7 64631 


F 


DEC 




.3 


F equals motor friction 


46 3 14631 












17 7 56314 


G 


DEC 




.2 


G equals motor inertia 


63 1 46314 












20 0 06314 


H 


DEC 




.8 


H equals load inertia 


63 1 46314 












20 0 14000 


ONE 


DEC 




1 . 


Integer 


00 0 00000 












00 0 00000 


GU 


DEC 




0 


Constant initially zero 


00 0 00000 












17 7 64631 


BACKLASH 


DEC 




.3 


Amount of backlash 


46 3 14631 












01 0 00000 


FACTOR2 


OCT 


0100000000000000 


00 0 00000 












00 0 00000 


INDEX 


OCT 




0 




00 0 00000 












00 0 00000 


INCRONE 


OCT 




1 


Used in count 


00 0 00001 












00 0 00000 


BLANK 


OCT 




0 




00 0 00000 












00 0 00000 


TEN 


OCT 




10 


Used in count 


00 0 00010 












17 7 15075 


1 


DEC 




.01 


Beginning of floating- 


34 1 21727 










point conversion table 
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20175 

20176 

20177 

20200 

20201 

20202 

20203 

20204 

20205 

20206 

20207 

20210 

20211 

20212 

20213 

20214 

20215 

20216 

20217 

20220 

20221 

20222 



TABLE I - Continued 



17 7 
34 1 


25075 

21727 


2 


DEC 


.02 


17 7 
12 1 


27534 

72702 


3 


DEC 


.03 


17 7 
34 1 


35075 

21727 


4 


DEC 


.04 


17 7 
63 1 


36314 

46314 


5 


DEC 


.05 


17 7 
12 1 


37534 

72702 


6 


DEC 


.06 


17 7 
60 5 


44365 

07534 


7 


DEC 


.07 


17 7 
34 1 


45075 

21727 


8 


DEC 


.08 


17 7 
07 5 


45605 

34121 


9 


DEC 


.09 


17 7 
63 1 


46314 

46314 


10 


DEC 


.10 


17 7 

36 5 


47024 

60507 


11 


DEC 


.11 


17 7 
12 1 


47534 

72702 


12 


DEC 


.12 


17 7 
72 7 


54121 

02436 


13 


DEC 


.13 


17 7 
60 5 


54365 

07534 


14 


DEC 


.14 


17 7 

46 3 


54631 

14631 


15 


DEC 


.15 


17 7 
34 1 


55075 

21727 


16 


DECq 


.16 


17 7 
21 7 


55341 

27024 


17 


DEC 


.17 


17 7 
07 5 


55605 

34121 


18 


DEC 


.18 


17 7 

75 3 


56050 

41217 


19 


DEC 


.19 


17 7 
63 1 


56314 

46314 


20 


DEC 


.20 


17 7 
36 5 


57024 

60507 


22 


DEC 


.22 


17 7 
12 1 


57534 

72702 


24 


DEC 


.24 


17 7 
00 0 


64000 

00000 


25 


DEC 


.25 
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20223 

20224 

20225 

20226 

20227 

20230 

20231 

20232 

20233 

20234 

20235 

20236 

20237 

20240 

20241 

20242 

20243 

20244 

20245 

20246 

20247 



TABLE I - Continued 



17 


7 


64121 


26 


DEC 


.26 


72 


7 


02436 








17 


7 


64365 


28 


DEC 


.28 


60 


5 


07534 








17 


7 


64631 


30 


DEC 


.30 


46 


3 


14631 








17 


7 


65075 


32 


DEC 


.32 


34 


1 


21727 








17 


7 


65463 


35 


DEC 


.35 


14 


6 


31463 








17 


7 


65605 


36 


DEC 


.36 


07 


5 


34121 








17 


7 


66314 


40 


DEC 


.40 


63 


1 


46314 








17 


7 


67024 


44 


DEC 


.44 


36 


5 


60507 








17 


7 


67146 


45 


DEC 


.45 


31 


4 


63146 








17 


7 


67534 


48 


DEC 


.48 


12 


1 


72702 








20 


0 


04000 


50 


DEC 


.50 


00 


0 


00000 








20 


0 


04314 


55 


DEC 


.55 


63 


1 


46314 








20 


0 


04631 


60 


DEC 


.60 


46 


3 


14631 








20 


0 


05146 


65 


DEC 


.65 


31 


4 


63146 








20 


0 


05463 


70 


DEC 


.70 


14 


6 


31463 








20 


0 


06314 


80 


DEC 


.80 


63 


1 


46314 








20 


0 


07146 


90 


DEC 


.90 


31 


4 


63146 








20 


0 


14000 


100 


DEC 


1.00 


00 


0 


00000 








20 


0 


14631 


120 


DEC 


1.20 


46 


3 


14631 








20 


0 


15463 


140 


DEC 


1.40 


14 


6 


31463 








20 


0 


16314 


160 


DEC 


1.60 


63 


1 


46314 
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TABLE I - Continued 



20250 


20 0 17146 


180 


DEC 


1.80 




31 4 63146 








20251 


20 0 24000 


200 


DEC 


2.00 




00 0 00000 








20252 


00 0 00000 


BUF 


BSS 


6 




00 0 00000 












DECO 


EQU 


71000 






RUNGE 


EQU 


60200 



END 
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1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 



TABLE II 



TYPICAL PRINTOUT OF DATA 



MOTOR 

SPEED 

77040630 - 1 
15278417 
22518046 
29387916 
35856424 
41895488 
47480616 
52590944 
57209252 
61321964 
64919110 
67994283 ’ 
70544563 
72570424 
74075629 
75067101 
75554775 
75551450 
75072608 
73561932 
70007602 
64634564 
57670572 
49317406 
39753812 
29138106 
17610490 
52951014 - 1 
- 76981710 - 1 
-21272070 
48719261 
44993108 
39850916 
34619795 
29335478 
24032905 
18746033 
13507650 



MOTOR 

POSITION 

30333708 - 2 
12247439 - 1 
27389548 - 1 
48177680 - 1 
74303153 - 1 
10543340 
14121485 
18127574 
22522912 
27267569 
32320675 
37640704 
43185760 
48913847 
54783141 
60752245 
66780440 
72827914 
78855989 
84814975 
90570699 
95967717 
10086982 1 

10515798 1 

10872835 1 

11149055 1 

11336614 1 

11428723 1 

11419529 1 

11304005 1 

14276023 1 

14557350 1 

14896797 1 

15194727 1 

15450572 1 

15664046 1 

15835140 1 

15964112 1 



LOAD 

SPEED 

68901844 - 1 

14489345 

21757968 

28660703 

35165619 

41244294 

46871885 

52027169 

56692568 

60854140 

64501553 

67628038 

70230320 

72308523 

73866071 

74909556 

75448598 

75495690 

75066030 

74326600 

73587038 

72854835 

72129917 

71412213 

70701649 

69998156 

69301663 

68612100 

67929398 

67253489 

48719261 

45474833 

40349320 

35131214 

29856282 

24559521 

19274963 

14035491 



LOAD 

POSITION 

30007982 - 2 
12215856 - 1 
27359124 - 1 
48148568 - 1 
74275496 - 1 
10540733 
14119047 
18125317 
22520843 
27265695 
32319002 
37639237 
43184500 
48912797 
54782300 
60751613 
66780013 
72827688 
78855961 
84831950 
90748447 
96606073 
10240541 1 

10814705 1 

11383155 1 

11945950 1 

12503145 1 

13054795 1 

13600957 1 

14141684 1 

14276023 1 

14557542 1 

14896996 1 

15194932 1 

15450780 1 

15664257 1 

15835352 1 

15964323 1 
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APPENDIX A 



DEVELOPMENT OF CHART PARAMETERS 



The purpose of this appendix is to present a procedure whereby the 
factors required to enter the design charts of this thesis may be computed 
for a given servo system. Equations for the field control motor, the arma- 
ture control motor, and the two phase motor are developed. The application 
of a given system represented by a differential equation of the form 



since the design charts are based on an undamped natural frequency of one. 
This can be accomplished by a time scaling of the problem. Thus, in the 



When these substitutions are made there results the desired form suitable 
for use with the design charts. From this , it is evident that the damping 



charts , while based on an undamped natural frequency of one, are applicable 
to other frequencies as well. 

The second order servo can be represented by a block diagram as 
follows . 




requires a transformation to the form 



0 



0 



first equation above if the independent variables are related by ' 

where 




and 



coefficient, , remains unaltered in the transformation, and the design 
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Three types of motors are considered, the armature control motor, the 
field control motor, and the two phase motor. For the period of time when 
the gears are in contact the system is a simple second order system. The 
differential equations for this combined motion can be derived by equating 
the load torque required to the motor torque developed. For the armature 
control motor this results in the following equation. 















For the field control motor if the field time constant is negligible the 



equation is similar. 

In like manner the differential equation for a two phase motor is 



ZZ y~ZZ 

where 



Ki = motor torque constant 

K , = motor generator constant 

= error measurement constant 
R = armature resistance 
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f* = gear constant 

T , = inertia of motor referred to load shaft 

M 

J = inertia of motor 

J m 

J = inertia of load 

L 

f m = friction constant of motor 

= friction constant of motor referred to load shaft 
F = friction constant of load 

Fj = friction of load plus motor friction referred to load shaft 
v _ £3 



(from characteristic curves of two phase motor) 
(from characteristic curves of two phase motor) 



r/? 



K 



N 



9T 

Bn 



K, 



ar 
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Procedure for Armature Control Motor 



1. 

2 . 

3. 

4. 

5. 

6 . 

7. 

8 . 

9. 



Find K from 



/<r= 



X/ X-7 



Find J referred to load shaft . 
M 



Find uJ, 



77 



^ Ik l 



Juf-S-J/L. 



Kn, 



Find F t from /T_ =- ^ y. 



Find F^ from 
Find 4 
Find J L /J M 




Find F L A Total = 







Use 6,7, and 8 to select proper design chart and find size of limit 
cycle for 0.3 rad backlash. 



10. Find size of limit cycle for desired backlash from 



/ .) kSkkkzkk^- ) 

A/a j/7~ ^ a// A/a//, / 



Procedure for Field Control Motor 



1. through 4. Same as Armature Control Motor 

5 . Find 4 ^ KOJn) 

6. Continue as Armature Control Motor 



Procedure for Two Phase Motor 



1. Find K from characteristic curves of motor 
e 



2 . Find T , referred to load chart 

M 

3 . Find Oj n = —Kk__ 



4. Find F 



total 



5. Find K from characteristic curves of motor, 
n 
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7. Continue as Armature Control Motor. 
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APPENDIX B 



SAMPLE PROBLEM 

A Elinco FD-162 115 volt D. C. Shunt Motor, 1/125 Horsepower 

4,000 Rpm is coupled through a gearbox to a load. The constants of the 

proportional error servo are as follows: 

Armature resistance, R — 282 ohms 

Motor Generator constant, K^O.233 volts/rad/sec 

Torque Constant, K«28.2 ounce - inch/ampere 

2 

Motor Inertia, J == 0.61 ounce -inch 
m 

-s 

Motor Friction, f —.00064 ounce-inch/rad/sec 
m 

Error measurement constant, Kg=10 volts/rad 
Gearbox factor, -^>=2 
Backlash 0 . 1 rad 

-6 2 

Load Inertia J =16.4 x 10 slug ft 

L 

Load Friction Fj^.002 ounce-inch/rad/sec 



Solution : Using the procedure for the armature control motor in appendix A. 
1 //= ^ /r-r- 



2 . = -4?r = 



0m £> / 






— 3. 2 x/o C Jx-c/G 



all — \ 1 / & £&d/sec 
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6 . 



> /-7b ±JlL V-L- j 1/^ __ 

> 5 ^^/ 6&.6*,oyc<zX'4*J - /9 



/_ 

r%/,} 



/70 x/O 



7. 



J2" __ 7.7 x//?' 

77 S.7 y /&'<=> 






= 7. a 



8 . 



^2 

/Z T577^ 




— .077' 



9 . From Figure 7 representing a friction ratio of . 1 , and for an 
inertia ration of 2.0 and a damping coefficient of .19 the 
size of the limit cycle for 0.3 radian backlash is .25 radian. 

10. The size of the limit cycle for 0.1 radian backlash is 
limit cycle = 07 — - *£3 . 

0 3/ - ■ 
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APPENDIX C 



USE OF DESIGN CHARTS 

Figures 7 through 11 are presented as representing design charts 
for the second order servo with backlash. These charts can be used to 
rapidly determine the size of the limit cycle by the following procedure . 

1. If the system variables are not known, compute these by 

the methods of Appendix A. 

2. Choose the chart representing the Load Friction/Total Friction 

desired. 

3. Read the Size of the Limit Cycle resulting from the intersection 

of the abscissa representing inertia ratio and the desired 
damping coefficient, interpolating if necessary. 

4. The value determined represents that Limit Cycle resulting 

from a backlash of 0.3 radians. For a given backlash the 
limit cycle is directly proportional, thus divide the Limit 
Cycle found in 3 above by 0.3 and multiply by the desired 
backlash. 
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